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Generating Detector Efficiency Estimates for a PET Scanner 

Field of the Invention 

This invention relates to a method of, and computer software for, 
5 generating detector efficiency estimates for a positron emission tomography 
(PET) scaimer including a detector array for generating detection data and a 
single photon source. 

Background of the Invention 

10 There are various different types of acquisition procedure in PET, 

including an emission scan^ a transmission scan and a blank scan. 

A typical emission scan begins with the injection of a solution including 
a tracer, which is a pharmaceutical compound including a radio-isotope with a 
short half-life, into the subject. The subject may be human or animal. The tracer 

15 moves to, and is typically taken up, in one or more organs in the subject 
according to biological and biochemical processes which occur within the 
subject. When the radio-isotope decays, it emits a positron, which travels a short 
distance before annihilating with an electron. This annihilation produces two 
high energy photons propagating in opposite directions. The PET scanner 

20 includes a photon detector array arranged (usually in a ring-shaped pattern) 
around the scanning area. If two photons are detected within a short timing 
window, a so-called "coincidence" is recorded along a line of response (LOR) 
connecting the two detectors. Coincidence counts along each LOR are 
incremented in a data storage part referred to as a sinogram bin each time a 

25 coincidence is detected between the corresponding detector pair. The output 
coincidence coxmt data firom each sinogram bin is typically processed using 
tracer uptake models and image processing techniques to obtain volumetric 
medical images and volumetric tracer uptake rate data for the subject. 

It would be desirable to be able to ascertain detector efficiencies in order 

30 to be able to detect deterioration in the state of operation of detectors in the 
detector array, since upon deterioration of a detector, the scan results are 
accordingly reduced in accuracy, leading to poor image quality. 
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For quantitative results from PET images, attenuation correction forms 
one of the most crucial data correction stages. In a conventional PET scaimer, 
the scanner is provided with one or more positron emitter rod sources, formed of 
a material such as ^^Ge, which emit dual aimihilation photons. To derive 

5 attenuation factors, two acquisitions using the rod sources are conventionally 
used — a blank scan, in which tiie subject being scanned is not present in the 
scanning area (typically, the scaimer is empty except for the presence of the 
sources) and a transmission scan in which the subject is present in the scanning 
area. Since the source material is a positron emitter, the two photons arising 

10 from the aimihilation of a positron and an electron are acquired in coincidence, in 
the same maimer as with an emission scan. Conventionally, the results of the 
blank scan are then divided by the results of the transmission scan to derive an 
attenuation sinogram. The attenuation sinogram is then used to correct the 
emission scan for attenuation, although extra processing steps on blank scan data 

15 and transmission scan data are sometimes used. 

Also, by collecting the coincidence count data during a blank scan and 
processing the count data using an efficiency estimation algorithm, individual 
detector efficiencies can be estimated. For non-rotating PET scanners, estimates 
of the detector efficiency may be calculated by adding together the counts 

20 detected between a particular detector and all the detectors that are in 
coincidence with it. This is known as a fansum and it is based on the following 
measurement model: 

Mij = £i€jAij (1) 

where My are the measured coincidence coimts between detectors i andy, 
25 and €/ are the intrinsic efficiencies of detectors i and j respectively and Ay are 
the ideal coincidence counts measured if the detectors had an ideal performance. 
Thus the measured counts Mfj of a particular coincident pair of detectors are 
proportional to the product of the individual detector efficiencies. An estimate of 
each individual detector efficiency €/ is made by smnming the measured 
30 coincidence counts My over all detectors j in coincidence with detector i: 
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The above equation is based on the assximption that: 

^SjAy = constant (3) 
J 

The assumption is reasonable providing that the eflBciencies are not 

significantly different, that J] Ay is roughly independent of i and that the number 

J 

5 of coincident detectors is sufficiently large. More accurate methods have been 
proposed which take into account noise and other variables. Further examples of 
efficiency estimation algorithms for calculating detector efficiencies in 
conventional PET scanners are described in **Maximum-Likelihood Estimation 
of Normalisation Factors for PET", D. Hogg, K. Thielemans, T. Spinks and N. 
10 Spyrou, Nuclear Science Symposium and Medical Imaging Conference Record, 
ffiEE, 2001. 

In an alternative PET scanner arrangement, used for example in the 
ECAT EXACT3D™ PET scanner, a single photon source is provided for the 
transmission scaiL In this scanner the source, a ^^^Cs pellet, is automatically 

15 moved by a transport mechanism having a spiral tube that is placed in firont of 
the detectors when a transmission scan takes place. By forcmg a liquid through 
the tube the source is moved at a constant speed througfh the scanner in a helical 
motion to provide a full 3D transmission scan. 

As it uses a single photon source however, there is no second photon 

20 generated during an annihilation with which to form a coincident pair. With no 
coincident pair, no real coincidence coimts are available from which a 
transmission sinogram could be created. The way this problem is solved is that 
the scanner in this mode measures only the detections furthest from the source, 
ignoring those registered in detectors close to the source, while tracking the 

25 position of the source. The source position is measured at four equidistant points 
about the helix. Since the speed of movement of the source is effectively 
imiform, its position can be calculated accurately to within one detector position 
over the entire length of the helix. From the position of the detector and fhe 
source, an artificial coincidence count is generated. The set of artificial 

30 coincidence counts detected during a transmission scan are collected and a 
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transmission sinogram is calculated. However, detector efficiency estimates are 
not obtained from these residts. 

Some PET scaimers have detector arrangements which rotate during an 
acquisition. Typically, the detectors are arranged in two or more banks which do 

5 not fully surround the subject. Altematively, the detectors may be arranged in a 
non-ring-shaped pattern. In these cases, the model given by equation (1) above 
is no longer appropriate, as the rotation causes the direct inter-relation between 
detectors and entries in the sinogram bins to break down. Every sinogram bin 
typically contaiiis counts detected by multiple detector-pairs. For some of these 

10 scanners, including the ECAT ART™ PET scanner, the transmission 
measurement is again performed using a single photon point source. However, 
the source is transported only axially relative to the (rotating) detector arrays. 
The ECAT ART scanner contains two such sources, one located on each of the 
opposing detector banks. The ECAT ART scanner is described in further detail 

15 in the article 'The ECAT ART Scaimer for Positron Emission Tomography: "1. 
Improvements in Performance Characteristics", David W. Townsend et al.. 
Clinical Positron Imaging, Vol. 2, No. 1, 1999 and the article "Design and 
Performance of a Single Photon Transmission Measurement for the ECAT ART" 
C.C. Watson, W.F. Jones, T. Brun, K. Baker, K. Vaigneur and J. Young, IEEE 

20 Medical Imaging Conference Record M9-02; 1998. 

A phantom scan involves positioning a body, referred to as the phantom, 
containing a positron-emitting radioactive source material in the scatmer. A 
phantom scan is either an emission or a transmission scan with a ^phantom' in the 
scaimer instead of a patient. Typically, the "phantom' is an object made out of 

25 plexiglass or suchlike and filled with water mixed with a radioactive substance. 
The phantom has known shape and attenuation characteristics. A phantom 
emission scan has in the past been lised to calculate detector efficiency estimates. 
A problem with having to conduct a phantom scan is that the phantom must be 
handled by an operator to place the phantom inside the scanner and to 

30 subsequently remove the phantom. It is both inconvenient and time-consimMng, 
since in any event an operator needs to be available to handle the phantom at the 
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start and end of the procedure. It is also potentially hazardous, in terms of lifting 
the phantom and exposure to radioactivity. 

It would be desirable to be able to estimate detector efficiency data in 
PET scanners of a type which use a single photon source for transndssion scans, 
5 without the drawbacks of using a phantom scan. 

In accordance with one aspect of the present invention there is provided a 
method of generating detector efficiency data for a positron emission 
tomography scanner including: 

a detector array for generating detection data; and 
10 a single photon source, 

wherein the method comprises: 

conducting an acquisition procedure using the single photon source to 
produce detection data; and 

processing said detection data using an efficiency estimation algorithm to 
15 calculate data representative of the efficiencies of individual detectors in said 
array. 

By enabling the estimation of detector efficiencies from an acquisition 
using the single photon source, detector efficiencies can be generated without 
significant inconvenience to an operator. The detector efficiencies may be 
20 derived firom a blank scan acquisition conducted at the operator's convenience. 
Furthermore, the regular need for the use of a phantom scan procedure can be 
avoided. 

Jn a preferred embodiment of the invention, detector efficiency estimates 
are made using the artificial coincidence coimts generated during a blank scan 

25 acquisition made using the single photon source. Typically, artificial 
coincidence counts are tiie only suitable detection data made available as an 
output fix>m a scanner of a type such as an ECAT EXACT3D PET scanner. 
However, detector efficiency data cannot be accurately estimated firom artificial 
coincidence counts using known techniques, because the known measurement 

30 models do not apply. The present invention provides a new measurement model 
and exemplary efficiency estimation algorithms, which can be applied to 
artificial coincidence counts produced using blank scans. 
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Further features and advantages of the present invention will become 
apparent from the following description of preferred embodiments of the present 
invention, made by way of example only with reference to the accompanying 
drawings. 

5 

Brief Description of Drawings 

Figure 1 is a schematic diagram of a PET scaiming facility; 

Figure 2 is a schematic diagram of data processing components in a non- 
rotating PBT scaxmer and an associated data processing terminal; 
10 Figure 3 is a schematic cross section of the arrangement within a non- 

rotating PET scanner during an emission scan acquisition; 

Figure 4 is a schematic cross section of the arrangement within a non- 
rotating PET scanner during a blank scan acquisition using a single photon 
source; 

15 Figure 5 is an illustration of a user interface provided in an embodiment 

of the invention for presenting detector efficiency estimates to an operator; and 

Figure 6 is a schematic cross section of the arrangement within a rotating 
PET scanner during a blank scan acquisition using a single photon source. 

20 Detailed Description of Preferr ed HmhnHinrientfi 

Referring now to Figure 1, a PET scaiming facility, arranged in 

accordance y/ifh an embodiment of the invmtion, is shown, which includes a 

PET scanner 2, a tracer generator module 4 and an operator computer terminal 6. 

The scanner 2 includes a detector array 8 arranged about a scanning area, in 
25 which a subject 10 is located during a transmission scan and during an emission 

scan. 

As shown in Figure 2, the PET scanner 2 includes a control imit 30, 
detection data processing circuitry 32, one or more single photon source position 
detectors 34, a count store 36 and an Input/Output (I/O) device 38. The 
30 computer terminal 6 includes a central processing xmit (CPU) 42, memory 44, 
hard disc drive 46 and I/O device 40, which facilitates intercoimection of the 
computer 6 with the PET scaimer 2. Operating system programs 48 are stored on 
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the hard disc drive 46, and control, in a known naanner, low level operation of 
the computer terminal 6. Program files and data 50 are also stored on the hard 
disc drive 46, and control, in a known manner, outputs to an operator via 
associated devices. The associated devices include a display, a pointing device 

5 and keyboard (not shown), which receive input fix>m, and ou^ut information to, 
tiie operator via further I/O devices (not shown). Included in the program files 
SO stored on the hard drive 46 are a detector efi&ciency calculating application 52 
and a detector efiBciency user interface application 54. A database 56 is used to 
store the detection data transferred firom the PET scanner 2. 

10 The detection data processing circuitry 32 of Ihie PET scanner op^tes in 

two different modes. In a positron source mode, the detection data processing 
circuitry 32 processes all events detected in the detector array 8, and by using a 
coincidence timing window, detects coincidences betwcCTi the events which are 
recorded as real coincidence counts Mij, These real coincidence counts are then 

15 output to sinogram bins in the count data store 36 for subsequent processing. 

In a single photon source mode, the detection data processing circuitry 32 
generates artificial coincidence counts when a single photon source is exposed in 
ibQ scanner, as will be described below in fiulher detail. The one or more single 
photon source position detectors 34 provide a position detection output to control 

20 unit 30 which generates a dynamic position estimate during tiie acquisition, 
which is sent to the detection data processing circuitry 32. The detection data 
processing circuitry 32 filters the detector outputs so as to discard all events 
detected on the side of the detector array nearest the current source position. 
Next, the detection data processing circuitry 32 generates an artificial 

25 coincidence count M'ij^ and transmits the artificial coincidence count data to 
siuogram bins in the coimt data store 36 for subsequent processing. 

One implementation of the invention, which is described in fiirther detail 
below, relates to a non-rotating PET scanner. In the embodiment of this 
implementation, a PET scanner of the ECAT EXACT3D™ type is used. 

30 However, it should be noted that this implementation of the invention applies 
also to other non-rotating scanners which include a single photon source. In this 
embodiment, the detectors in the detector array 8 are arranged in square detector 
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blocks, each containing 64 detector elements. The detectors are arranged in six 
rings of detector blocks. There are 48 detector blocks in each ring. The total 
nxunber of lines of response (LORs) that can be acquired is thus, in this 
embodiment, about 190 milUon. 
5 Figure 3 illustrates features of operation of the non-rotating PET scanner 

in a positron source mode during an emission scan. During an emission scan, the 
subject 10 is placed in the scanning area and contains the tracer generated by the 
tracer g^erator module 4. Figure 3 shows a positron emission event being 
registered in the detector array 8. The positron 60 annihilates and generates a 

10 first photon 62 travelling in one direction and a second photon 64 travelling in an 
opposite direction. The first photon is detected by a detector element 70 on one 
side of the detector array 8, and the second photon 64 is detected in a different 
detector element 72 on the other side of the detector array 8. As described 
above, the detection data processing circuitry 32 registers the two as a 

15 coincidence along the LOR defined between the two different detectors 70, 72. 

Figure 4 illustrates features of operation of the non-rotating PET scanner 
in a single photon source mode during a blank scan. In this mode, a single 
photon source 80 is transported around the periphery of the scanning area using a 
single photon source transport mechanism 82. The single photon source 80, 

20 referred to as a **poinf ' source due to its small diameter, consisting of a pellet of 
^^^Cs is driven in a fluid-filled steel tube wound into a helix and positioned just 
inside the detector array 8. The tube diameter is 4mm or smaller and the *^^Cs 
point source, with a radioactivity typically between 1 to 10 mCi, has a diameter 
of 2mm or less. The source is transported in a Hquid carrier and is buoyancy 

25 matched to the specific gravity of the liquid. The arrangement is designed to 
mcrease consistency of the speed of travel of the source. 

There are 48 turns of the helix within the detector field of view, 
corresponding to one turn per detector ring. The source carrier is driven by a 
pump (not shown) at a speed of about Ims"^ The source 80 is shown in Figure 4 

30 emitting a single photon 90 which is detected at a detector 94 on the opposite 
side of the detector array 8 to the source 80. Figure 4 shows a virtual second 
photon 92, representing the artificial coincidence registered in the detector data 
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processing circuitry 32 as described above. An event along an LOR is registered 
between the two respective detectors involved, the detector 94 registering the 
event and the detector 96 immediately adjacent the current source position 80. 
As described above, the detection data processing circuitry 32 generates an 
5 artificial coincidence count which is transmitted to the count data store 36. 

After the end of the acquisition, the operator can instruct computer 
terminal 6 to download the detection data held in count data store 36 to the 
database 56 on the computer terminal hard drive 46, where the data can be 
processed. 

10 In accordance with this embodiment of the present invention a 

measurement model for an acquisition using the single photon source is 
provided, which takes into account the method whereby the artificial coincidence 
counts are produced. The measurement model has the form: 

M^y^eiA\j^sjA^ji (4) 

15 where M'y are the measured artificial coincidence counts between detectors i and 
7, €/ and £j are the intrinsic efficiencies of detectors i and j respectively and A *ij 
are the ideal artificial coincidence counts measured if tibe detector efficiency of 
detector i had an ideal performance, while the detection efficiency of detector j 
was 0. This measurement model takes into account that an event M'ij recorded 

20 between detectors i and J will have originated either firom detector i or j\ These 
will have each have a respective inherent efficiency e,- or €j. For a ring-shaped 
scanner where the single photon source orbits a spiral with a centre of rotation 
coinciding with the axis of the PET seamier, it can be assiraied that A A 'ji. 

Thus, the measurement model used is additive, in that the measured 

25 artificial coincidence counts M'y of a particular artificially coincident pair of 
detectors are equal to a weighted sum of their individual efficiencies. This 
contrasts with the known model for a tme coincidence system, based on equation 
(2) above, A^s^ch is multiplicative, in that the measured counts My of a particular 
coincident pair of detectors are proportional to the product of the individual 

30 detector efficiencies, which can lead to very significantly different results. 
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Next, an algorithm is provided for estimating detector efficiencies. 
Assuming that the transmission scan is reasonably constant over the region of the 
coincident pairs taken into consideration, then the A terms can be ignored. 
Equation (4) can then be rearranged and solved as an efficiency estimation 
5 algorithm using the following iterative formula: 

=T|^Z^V-Z^i] Niter (5) 

where J is the number of detectors in the accepted fan and k is the iteration 
number. For the first iteration the efficiencies can be initialized to any positive 
constant. For a fast convergence to a solution, new estimates of are passed 

10 into the right hand side of equation (5) as soon as they have been calculated. If 
the A 'ij terms cannot be ignored, a general iterative linear equation solver may be 
used. It is possible to use more sophisticated methods, for example based on a 
Maximum Likelihood prescription which takes noise in the measurements into 
accoxmt. In particular, because of the additive nature of the measurement model 

15 (4), the Maximum LikeUhood methods applicable to image reconstruction in 
PET, such as the MLEM algorithm as described in the article '^Maximum 
likeUhood reconstruction for emission tomography", Y.Vardi, L.A Shepp, IEEE 
Trans, Med. Imag., Vol. 1, pp.113-121,1982, may be used. To briefly describe 
the MLEM algorithm, we introduce some notation as follows. We use the index 

20 z to run over all sinograms bins. We define a matrix P to concisely write the 
linear measurement model as: 

^z-llP^^j (6) 
J 

or in matrix notation: 

M=P€ (7) 

25 The MLEM algorithm can then be written as: 

=(5* ^Pu)^P'\f^^Pe^) (8) 
where u is a sinogram with all elements equal to 1, -^ denotes element-wise 
division and 0 denotes element-wise multiplication. The initial estimate for s is 
arbitrary, as long as all elements are strictly positive. However, we recommend 
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setting all sj to 1. For more details, we make reference to the above mentioned 

article, the contents of which are incorporated herein by reference. 

Whichever efficiency calculating algorithm is selected, the detector 
efficiency calculating component 52 implements the efficiency calculating 

5 algorithm as provided hereby. 

A PET scanner can perform a transmission scan in either 2D or 3D mode. 
For patient scans on the ECAT EXACT 3D, both 3D blank and 3D transmission 
acquisitions can be conducted. Data from several LORs are then combined 
together in a process called panning to save data storage space. This means that 

10 the above measurement model (4) is no longer appropriate. To avoid this 
problem, a blank scan acqxiisition is preferably conducted in a 2D mode (in 
which no spanning occurs), and from this individual detector efficiencies can be 
calculated. However, where it is possible to acquire the blank scan acquisition 
data in 3D mode without spanning the data, the 3D blank scan data may be used 

15 to calculate detector efficiencies in accordance with the invention. Even with 3D 
blank scan data witii spanning, although more complicated, the measurement 
model is still a linear combination of the detector efficiencies, and hence this 
system of linear equations can be solved by an iterative procedure similar to 
equation (S) above. Similarly, a Maximum Likelihood approach such as the 

20 MLEM algorithm can be sailed. 

To experimentally test the efficiency estimate algorithm provided using 
equation (5) above, a direct plane 1 hour blank scan acquisition was performed 
on an ECAT EXACT3D scanner. The efficiencies were calculated using 
equation (5) and a corrective efficiency sinogram generated for each ring. The 

25 acquired blank sinogram, generated using the artificial coincidence coimt data, 
showed clearly the differing detector efficiencies by diagonal lines of varying 
intensity. The acquired blank sinogram was corrected for these effects using the 
calculated efficiency sinogram, and the corrected blank scan sinogram was seen 
to be free of diagonal lines of varying intensity. Thus, the calculated efficiency 

30 estimates were shown to be accurate. 
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Figure 5 illustrates a display part of a user interface generated on the 
computer terminal 6 by the detector eflBciency user interface apphcation 54 using 
the detector efficiency data calculated by data efficiency calculating application 
52 after a blank scan procedure as described above. The display includes a 
5 detector representation 100, showing each of the detectors in the detector array 8 
on a flat grid 100. The flat grid 100 is divided into detector block squares 102, 
representing each detector within the block as a separate pixel for the different 
brightness intensity. The brightness intensity used in each pixel represents the 
calculated efficiency estimate for the respective detector. A relatively dark pixel 

10 104 is used to illustrate a detector having a relatively poor eflBciency. A 
relatively bright pixel 106 is used to illustrate a detector having a good, fully 
operative efficiency. An intermediate brightness intensity pixel 108 is used to 
represent a detector having an intermediate efficiency. Thixs, by viewing the grid 
100, the operator is able to visually identify detectors having poor and 

15 intermediate efficiencies, thereby to detect problems in the detector array 8. For 
example, occasionally an entire detector block will be somehow impaired or 
even inoperative. In this case, the entire detector block square 102 will appear 
relatively dark, and the operator can determine that the block requires 
replacement in order to keep the PET scanner image quality high. 

20 Furthermore, the detector efficiency user interface application 54 is 

provided with one or more automated block impairment check routines which 
are used to highlight blocks which are determined to be below a given threshold 
of acceptable operation. In a first check routine, an average efficiency is taken 
across each block. If the average detector efficiency within any block is below a 

25 preset threshold, flie relevant block square 102 is highlighted, for example by the 
addition of a coloured border or a flashing display of the block square 102. 

In a second block check routine, the coefficient of variation of 
efficiencies within a block is calculated. The detector efficiency user interface 
application 54 then highlights, for example by the addition of a coloured border 

30 or a flashing display, any detector square 102 representing a block having a 
coefficient of variation which is higher than a preset threshold. 
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Other implementations of the invention, which are described in further 
detail below, relate to a rotating PET scanner. In an embodiment of this 
implementation, a scanner of the ECAT ART™ type is used. However, it should 
be noted that this implementation of the invention applies also to other rotating 

5 scanners which include a single photon source. In this particular embodiment, 
the detectors in the detector array are arranged in square detector blocks, each 
containing 64 detector elements. The detectors are arranged in two banks of 
detector blocks. The banks do not entirely surround the scanning area, thus 
allowing fewer detector elements to be used in the scanner. 

10 Figure 6 illustrates features of operation of the rotating PET scanner in a 

single photon source mode during a blank scan. The scaim^ has first and second 
detector banks 120, 122 arranged on difTerent sides of the patient The two 
detector banks are arranged on a rotating mechanism which allows the banks to 
rotate about the scanning area. For blank scans, a single photon source 124, 130 

15 is held in a single photon source transport mechanism 126, 132 on one side of 
each of the respective detector banks 120, 122. The single photon sources 124, 
130 are thus mounted to, and rotate with, the detector array. Each single photon 
source transport mechanism is arranged to transport the source axially of the 
rotating detector array during a blank scan. A single photon source shield 128, 

20 134 is provided for each single photon source. During a blank scan, one or each 
of the single photon source shields 128, 134 may be retracted to expose the 
single photon source in the interior of the scanner. Figure 6 shows an 
arrangement in which the first single photon source 124, on the first detector 
bank 120, is e?q)osed by the retraction of the first single photon source shield 

25 128, and in which flie second single photon source 130, on the second detector 
bank 122, is exposed by the retraction of the second single photon source shield 
130. In this arrangement, a photon 138 emitted by the first single photon source 
124 is detected at a detector element 140 on the second detector bank 122. 

For a rotating scanner, a count in a given sinogram bin corresponding to a 

30 single LOR, may correspond to different detector pairs, because when the 
detector array is rotating, the LOR between two detectors at time 0 would later 
be covered by another detector pair. 
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In one embodiment of the invention, detector efficiency estimates are 
obtained by conducting a blank scan acquisition when the rotatable detector 
array is held stationery. In the case that no LORs are combined (i.e. no 
spanning), there is a one-to-many correspondence between the (non-zero) bins in 
5 the sinogram and the detectors. In this case, the measurement model is: 

M\=ej^A\ (9) 
where there is a one-to-many correspondence between jz and and 
efficiencies can be solved by using: 

10 where the set Zj contains all sinogram .bins that correspond to the 

detector^', i.e. 

Zj^{z\J,^j} (11) 
Equation (10) is the Maximum Likelihood solution for the efficiencies 
when the measurement model (9) applies. The sum in (10) over all sinogram 
15 bins that correspond to the same detector reduces the noise in the estimated 
efficiencies. 

In fiirther embodiments of implementation of the invention for rotating 
scanners, detector efficiency estimates are generated when the detector array is 

20 rotating. The detector efficiencies may depend on temperature, and hence 
rotation speed, thus detector efficiency estimates may be generated which more 
accurately reflect the operative efficiencies. 

In one embodiment in which the detector array is rotated during a blank 
scan acquisition, one single photon source is exposed to the (rotating) detector 

25 array, and the other single photon source is shielded, as shown in Figure 6. In 
this case, a measurement model that can be used has the same form as equation 
(9) above, albeit with different A\ values, where there again is a one-to-many 
correspondence between jz and z. Hence a prescription similar to (10) can also 
be used. The measurement is then to be repeated with the other source viewable 

30 such that the other detector bank is irradiated. 
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In a further embodiment in which the detector array is rotated during a 
blank scan acquisition, both single photon sources are exposed in the (rotating) 
detector array. In this embodiment, at least some of the sinogram bin values will 
correspond to either two or four overlapping but dissent source-detector 

5 arrangements. A possible solution to this problem is to use gating hardware. 
Gating hardware provides multiple sinograms (one for each gate), where counts 
are assigned to a gate according to the time they are detected. The overlap effect 
occurs at the periphery of the detector banks. There is an angle, 6, which defines 
the smallest rotation angle before ovCTlap occurs. Taking the largest angle 6* <& 

10 such that 360^ is an integer multiple of 6', i.e. 360*^ = 6' * iV, where AT is an 
integer, N gates are set up, each of length T/N where Tis the period of rotation of 
the scanner. After a full rotation, the scanner hardware will add counts to the 
same set of gates again, such that N different sinograms, each witii uniquely 
identifiable correspondence between the LOR and detector, are produced. A 

15 measurement model, and analysis function, such as that described above using a 
single exposed source may then be employed to generated detector efficiencies. 

Finally, also for the case of rotating scanners, spanning is often 
introduced to reduce storage demands. As for non-rotating scanners, this will 
mean that the measurement model (9) is to be modified such that the measured 

20 artificial coimts are a linear combination of efficiencies. As before, one can solve 
for the efficiencies with a Maximum LikeUhood algorithm such as the MLEM 
algorithm. 

Thus, it will be appreciated that the present invention may be 
25 implemented to generate detector efficiencies firom a blank scan that could then 
be used as a daily check on the performance of a PET scaimer. The invention 
may also be used to provide detector efficiency factors for a scanner 
normahsation procedure. Since detector efficiencies change over time, it is 
desirable to be able to ascertain detector efficiency factors on a regular basis and 
30 without the need for complex procedures. The above-mentioned efficiency 
estimation algorithms, or alternatives thereto, may be used in such a procedure. 
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The above embodiments are to be understood as illustrative examples of 
the invention. It is to be understood that any feature described in relation to any 
one embodiment may be used alone, or in combination with other features 
described, and may also be used in combination with one or more features of any 
5 other of the embodiments, or any combination of any other of the embodiments. 
Furthermore, equivalents and modifications not described above may also be 
employed without departing fi-om the scope of the invention, which is defined in 
the accompanying claims. 



